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Abstract 

We prove functional limits theorems for the occupation time process of a system of particles 
moving independently in M. d according to a symmetric a-stable Levy process, and starting off from an 
inhomogeneous Poisson point measure with intensity measure n(dx) = (1 + \x\' y )~ 1 dx, 7 > 0, and other 
related measures. In contrast to the homogeneous case (7 = 0), the system is not in equilibrium and 
ultimately it vanishes, and there are more different types of occupation time limit processes depending 
on arrangements of the parameters 7, d and a. The case 7 < d < a leads to an extension of fractional 
Brownian motion. 
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1 Introduction 

Several authors have studied systems of particles moving independently in M. d according to a 
Markov process (usually a symmetric a-stable Levy process, < a < 2), and also systems having in 
addition a branching mechanism (e.g. 00013110111111311111311311011111131131^1 and references 
therein) . A typical assumption in the cited references is that the system starts off from a homogeneous 
Poisson point measure, i.e., with intensity the Lebesgue measure (denoted here by A). This assumption 
represents a strong technical simplification because in the special cases usually studied A is invariant 
for the semigroup of the motion, and this implies that the particle system without branching is in 
equilibrium, and for d > a a critical branching system converges towards equilibrium |17j . In this 
case the systems have been extensively studied. New situations appear if the initial condition is an 
inhomogenous Poisson point measure. 

In this paper we consider the system without branching, with symmetric a-stable Levy process 
for the particle motion, and initial inhomogeneous Poisson point measure with intensity measure \x of 
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the form 

dx 

uldx) = i— ;— , 7 > 0, 

and other more general related measures. In this case the system is not in equilibrium and ultimately 
it vanishes (see below). Therefore one should expect different types of results from those of the 
homogeneous case. Our purpose is to obtain functional limits for the rescaled occupation time process 
of the particle system in different cases. 

The particle system is described as follows. Given a Poisson point measure on M d with intensity 
measure fj,, particles evolve from its atoms, moving independently according to a symmetric a-stable 
Levy process (called standard a-stable process). Let N = (N t )t>o denote the empirical measure 
process of the system, i.e., 

tft = !%(*)> (i-i) 

i 

where {xi(t)}i are the positions of the particles at time t. Note that N-t converges in probability to 
the null measure as t — > oo (Appendix). Let Xt = (-Xr(i))t>o denote the normalized occupation time 
fluctuation process of the system, defined by 

-Tt 



1 ( Lt 

X T (t) = — (N s - EN s )ds, (1.2) 

-FT Jo 



where T is the time scaling and Ft is a norming. The problem is to find Fj> such that the process Xt 
converges in distribution as T — » oo (i.e., the time is accelerated), and to describe the limit process X 
in the cases where it exists. 

In the homogeneous case (corresponding to 7 = 0), the occupation time fluctuation limit process 
has three different forms, for d < a [S], d = a. and d > a [Sj. In the inhomogeneous case there are more 
results depending on the values of 7 relative to d and a when j<d: ^<d<a,j<d = a,^<a< 
d,j = a<d,^ = d<a,~f = d = a,a<~f<d. For "small" 7, i.e., 7 < a, the results are analogous 
to those of the homogeneous case, while for "large" 7, i.e., 7 > a, and this seems unexpected, they 
are of a different kind. The case 7 < d < a leads to a long range dependence, self-similar, centered 
Gaussian process £ with covariance 

"SAt 



[•SAt 

EUs= / u a ((t-u) b + (s-u) b )du, 
J 



where a = — 7/a 6 (—1,0), b = 1 — 1/ a G (0, 1/2], which is an extension of fractional Brownian motion 
with Hurst parameter H £ (1/2,3/4] (corresponding to 7 = 0; the process £ with maximal ranges for 
the values of the parameters a and b is discussed in j^j). Nevertheless, although the process £ depends 
on 7, its dependence exponent [Zj is independent of 7. The cases 7 = a < d and 7 = d = a give a new 
type of limits (with no counterpart in the homogeneous case), namely, centered, constant (and hence 
continuous) Gaussian processes on (0,oo), discontinuous at 0. 

For 7 > d the measure \i is finite, and the results are in sharp contrast to those for 7 < d. In 
this case we give the results for a finite measure fi in general, and for d < a they are akin to the 
famous limit theorem of Darling and Kac for the occupation time (without centering), and its 
generalization to path space by Bingham |2j. 

All the occupation time limit theorems are formulated in the context of 5 / (M' i )-valued processes, 
where 5'(IR d ) is the usual space of tempered distributions (dual of the space 5(M rf ) of smooth rapidly 
decreasing functions). In some cases the limit process is of the form A multiplied by a real valued 
process, but in others the limit is "truly" 5'(IR rf )-valued. In all cases where the particle motion is 
recurrent (d < a), the spatial structure of the limit process is A, independently of 7. 
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The methods of proof for the fluctuation limit theorems are analogous to those developed in |B] , 
with some new technical complexities because the measure fi is not invariant for the semigroup of 
the motion. On the other hand, there is a significant difference in the tightness proofs, as they 
require estimates for moments of arbitrary high order (whereas in Ej order 2 or 4 was enough). 
For the results of Darling-Kac type we proceed similarly as 2 . However, in our setting the uniform 
convergence condition (A) for that kind of results is not satisfied, and some additional work is needed. 

Convergence in distribution in the space of continuous functions C([0, r], for any r > 
is denoted by In some cases the interval [0,r] is replaced by [e,t],0 < e < r, because the limit 
process is discontinuous at t = 0. 

The duality between the spaces S'(M. d ) and iS(R d ) is denoted by (•, •). 

Generic constants are written C, Ci, C2, ■ ■ ., with possible dependencies in parenthesis. 

In section 2 we present the results. Section 3 begins with an explanation of the general method 
used for the proofs of the occupation time fluctuation limits, and then we prove most of the results. 
Some proofs that are similar to others are omitted, with some comments. 

The branching particle systems in the inhomogeneous case produce fewer results, but there are 
other kinds of difficulties related to extinction of the system. These results will be presented elsewhere. 



2 Results 

Let N and Xt be the processes defined in Ql.lj) and (|1.2|) . As stated in the Introduction, for 
simplicity most of our results are formulated for \x of the form 

^dx) = - d * , , 7 > 0. (2.1) 

^ ) i + | x |7' ' v ; 

Note that [x is finite for d < 7. More general measures \i will be considered later in this section. In 
the theorems below, if is a number depending on a, d, /j,, which may vary from case to case, and may 
be computed explicitly in each specific case. 

Different arrangements of a,^y,d yield different results, and we order them according to the rela- 
tionship between 7 and d. We start with 7 < d. 

Theorem 2.1 Let 7 < d < a (hence d = 1) and 

Ft = T l ~( d+1 ^ 2a . (2.2) 

Then Xt as T — > 00 , where £ is a real centered Gaussian process with covariance 




This theorem is a generalization of Theorem 2.1 in p . which corresponds to 7 = 0. 

Remark 2.2 The following properties of the process £ are easy to obtain. 

(a) For 7 = 0,£ is a fractional Brownian motion with Hurst parameter 1 — l/2a [5]. 

(b) £ is self-similar with index 1 — (1 +7) /2a. This is immediate from (2.3), but more generally, 
from the form of the fluctuation process given by (1.2), it follows that if Ft has the form T K f(T), 
where / is a function slowly varying at infinity and k > (as it is in our cases), then the limit process 
is self-similar with index k . 
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(c) £ is a long range dependence process where 

E(£ 8+r -&+T)(£r-&) = 0(T- 1 /") as T^oc, 

for < r < u , < s < t . Note that the dependence exponent [7j 1/a is independent of 7. 

(d) £ is not a Markov process and not a semimartingale. The non-semimartingale property can be 
proved by Lemma 2.1 in [I]. 

The next two theorems are generalizations of Theorem 2.1 in 6 (for 7 = 0). 

Theorem 2.3 Let 7 < d = a (= 1 or 2) and 

F T = (TlogT) 1 / 2 ^/ 20 . (2.4) 
Then Xx ^ KX/3 as T —* oo, where (3 is an inhomogeneous real Wiener process with covariance 

Efh/3. = {t/ : S) , ■ (2.5) 
1 — 7/a 

Theorem 2.4 Ze£ 7 < a < d and 

i?T = T a-7/«)/2_ 

XTien ^ Klf as T —>■ 00, where W is an S' (M. d ) -valued time inhomogeneous Wiener process with 
covariance functional 

E(W(t),<p 1 )(W(s),<p 2 ) = (* As ) 1 7/ ° / ^(x)G^ 2 (x)^, ^2 e S(M d ), (2.7) 

where G is the a-potential operator, i.e. 

Gip(x) = C a>d [ V(V L a dy, (2.8) 
Ju d \x - y\ d 01 

with C a4 = r(^)(2 Q vrlr(f j)- 1 . 

The analogy with the case 7 = breaks down for "large" 7, i.e., 7 > a. 
Theorem 2.5 Let 7 = a < d and 

F T = (logT) 1 / 2 . (2.9) 

Then X T =¥ KX in C([e, r], «S'(Ir)) as T -> 00 /or an?/ < e < t, where X is an <S'(R ) -valued 
Gaussian process constant in time on (0, oo),X(t) = X(l), and X(l) is centered with covariance 
functional 

E{X(l),i Pl ){X(l),i P2 }= [ i p 1 (x)Gi P2 (x)dx, i Pl ,i P2 eS(R d ), (2.10) 
where G is given by \2 . 8\) . 

Note that the limit process is discontinuous at t = since Xt(0) = 0. 

To complete the case 7 < d it remains to consider a < 7 < d. It turns out, however, that if a < d 
and a < 7, then the relationship between 7 and d is irrelevant. In this case the total occupation time 
is bounded, so it does not make sense to investigate the fluctuation process. More precisely, we have 
the following simple proposition. 
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Proposition 2.6 Let a < d and a < 7. Then 

POO 

E (N s ,ip)ds < 00, tp£S(R d ). 
Jo 

We now proceed to the critical case 7 = d. 
Theorem 2.7 Let 1 = d = 7 < a and 

F T = T l ~^ a {\ogT) 1 / 2 . (2.11) 

Then Xt ^ KX^ as T — > oo ; where £ is as in Theorem 2.1. 

The next case is "doubly critical" . 

Theorem 2.8 Let 7 = d = a (= 1 or 2) and 

F T = (logT) 3 / 2 . (2.12) 

T/ien X T =>■ if A?? in C([e,t],S'(R )) as T —>■ 00 for any < e < r, uAere 77 zs a rea/ Gaussian process 
constant in time on (0, 00), rjt = rji, and rji is standard normal. 

Here, as in Theorem 2.5, the limit process is discontinous at t = 0. 

So far we were assuming that fx is of the form (|2.1|) . It is rather clear that for 7 < o! we can take 
fi(dx) = \x\~^dx. Moreover, a careful analysis of the proofs shows that in this case fi can have a more 
general form, given in the following proposition. (For the case 7 = d, see the discussion after the proof 
of Proposition 2.9). 

Proposition 2.9 All the previous results for the case 7 < d remain true (with possibly different 
constants K ) for an intensity measure \i of the form 

H(dx) = u (dx) + Y^j^dx, (2.13) 

where v is a finite measure, and h is a nonnegative bounded function such that there exists a strictly 
positive limit 

lim / h(x)dx. (2.14) 



R^oo R d J\ X \< R 

It seems interesting and perhaps unexpected that it is not sufficient to assume that n{dx) = g(x)dx 
with 

C ] , < g(x) < ° 2 , . (2.15) 

1 + \x\-y ' ~ 1 + \x\i y ' 

We have the following counterexample. 

Example 2.10 Let 7 < d < a (d = 1) and let [i be of the form (|2.13|) with v = 0, and 

{1 for \x\ < 4, 
1 for (2k) 2k < \x\ < (2k + l) 2k+1 , 
2 for (2k + l) 2fc+1 < \x\ < (2(k + l)) 2 ( fc+1 ), 

k = 1,2,.... The limit (|2.14|) does not exist for this measure, whereas (|2.15j) obviously holds. The 
only nontrivial normalization (cf. Theorem 2.1) is that given by (|2.2[) . but we will explain later that 
the corresponding Xt does not converge as T — > 00. 
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There remains the case 7 > d. Here the situation changes dramatically and the results are of an 
entirely different nature. In particular, they do not depend on 7, but only on the fact that the measure 
a is finite. Therefore we will formulate our results for a general finite measure p. It turns out that 
the appropriate normalization is 

F T = T x ~ lta (2.16) 

if 1 = d < a, and 

F T = logT (2.17) 

if d = a. 

It is easy to see that in both cases f Q (N s , ip)ds converges to a finite limit as T — > 00, ip G 5(IR d ), 
hence there is no reason to consider fluctuation processes and it suffices to investigate the occupation 
process 

1 f Tt 

Y T (t) = — / N s ds. (2.18) 

FT Jo 

For a given a > 1, let L denote the local time process (at 0) of a standard real a-stable process. 
See, e.g., [Q for properties of L. In particular, L is a continuous increasing process, L(0) = 0. 
The relation between the processes Yt and L is given in the following theorem. 

Theorem 2.11 Let 1 = d < a and n be a finite measure on K. Let L\,L2, . . . be independent copies 

of L and let v be a Poisson random variable with parameter /jt(W) independent of L±, L2, Then for 

F T defined by lETA) . 

as T — > 00. 

This theorem is based on the following lemma, which is of interest by itself. 

Lemma 2.12 Let a,d and be as before, let £ be a real standard a-stable process, and denote by 
Zt its normalized occupation process, i.e., 



Then 



as T — > 00. 



1 f 

(Z T (t),<p) = — <p(( s )ds, V eS(R d ),t>0. (2.19) 

b T Jo 

Z T => KLX (2.20) 



This lemma is closely related to the famous Darling-Kac result . Their theorem was generalized 
by Bingham [2], who obtained the limit in path space for fixed positive (p with compact support. 
Fitzsimmons and Getoor |16j mention this limit for fixed general <p. We will present an outline of a 
proof of the lemma in the next section. 

It remains to consider the case d = a. 

Theorem 2.13 Let d = a (= 1 or 2) and u be a finite measure. Let pi,p2,... be i.i.d. standard 
exponential random variables and v a Poisson random variable with parameter fJ,(M.), independent of 
PI,P2, Then for Ft defined by \2.11\j , 



/ j 1 



in C([e,T],S'{R d )) as T -» 00 for any < e < r. 
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So the limit process is constant in time on (0, oo). 



3 Proofs 



3.1 General scheme 



We describe a general method used in the proofs of Theorems 2.1-2.9. Given < a < 2, let % denote 
the transition semigroup of the standard a-stable process C in R d , i.e. %<p =p t *<p, where p t is the 
transition density of £ . 

It is well known that, by the Poisson property, 

E{N t ,tp)= I %<p(x)v(dx), tp£S(R d ). (3.1) 

JM. d 

For a continuous 5 / (IR c( )-valued process X we define an ^'(IR^ 1 ) random variable X by 

{X,&)= [ (X(t),$(-,t))dt, $G5(R d+1 ). (3.2) 
Jo 

As explained in in order to prove Xt ^> X, where X is the limit process occuring in the specific 
theorem, it suffices to show that 

(X T , $) =>{X,$), $ e S(M d+1 ), (3.3) 

and that the family {(Xt, </?)}t>2 is tight in C([0, r],R) for any r > 0, for each ip G 5(IR d ). 

This scheme should be modified in an obvious way if we consider C([s,r],S'(R )). Without loss 
of generality we will always assume r = 1. 

Since the limits are Gaussian, in order to obtain (|3.3|) it suffices to show that 

lim Ee-& T & = # e -<*'*> (3.4) 

for any nonnegative <1? 6 5(IR d+1 ) (see, e.g., [Sj). 
Given such <3? we denote 



(3.5) 



By JL2J), (f3~T|) . (I3~2l . we have 

(X T ,$)= f (N u ,y T (;u))du- f [ %* T (-,u)fi(dx)du. (3.6) 
Jo Jo Jm d 



Hence, by the Poisson property, 

i-T 



Ee -(X T ,$) =exp |^ J ^T u ^ T (-,u)(x)[i(dx)du^ex.p^- J d V T (x,T)fi(dx)^, (3.7) 



where 

rt 



v T (x,t) = 1- £expj- J ^ T (x + (u,T -t + u)du\ , < t < T. (3.8) 
(Recall that £ is the standard a-stable process). 
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We know that (repeating the argument in for V = 0), by the Feynman-Kac formula, vt satisfies 
v T (x,t)= [ T t _ s (y T (;T-s)(l-v T (;s)))(x)ds. (3.9) 

J 

We will often use an immediate consequence of ()3.8[) and (|3.9|) : 

v T {x,t)< I T t _ s V T (;T-s)(x)ds. (3.10) 

J 

Putting (|3.9|) into (|3.7jl and then using H3.9j) once more we obtain 

Ee ~{X T ,Z) =e I(T)-II(T)^ (3n) 

where 

I(T)= f [ T T - s (^t(;T-s) [ S T s . u ^ T (;T-u)du)(x)ds^dx) (3.12) 

JR d JO V JO / 

and 

J/(T) = / [ T T T _ s (y T (-,T-s) [ S T S ^ T (-,T -u)v T (-,u))du] (x)d Sf i(dx). (3.13) 

JR d JO V JO / 

To prove (|3.4|) we will show that 

lim e'W = Ee~&>*\ (3.14) 

T — >oo 

and 

lim II (T) = 0. (3.15) 

T — >oo 

For simplicity we will prove (|3.14|) and (|3.15|) for <I> of the form 

$(x, t) = ip(x)ip(t), <p e S(R d ), i/i e S(R), <p,ip>0. (3.16) 
For such $ it will be convenient to denote 

tb(s)ds. V'rm = vl 

T 



X (t) = I ip(s)ds, XT(t)=x[ I): (3.17) 



then 

y T (x,t) = ^(x) X T(t). (3.18) 

Note that expressions I(T) and II(T) have more complicated forms than those corresponding to 
7 = •"). (i . since the measure is not invariant under % and it is infinite, so in particular the Fourier 
transform technique we have used before is not applicable. 

In order to prove tightness of {(Xt,<p)}t>2 for a given <p G 5(IR d ),c/? > (it suffices to take ip 
nonnegative) , we need a formula for the Laplace transform of (Xt^) — Xxiti), <p) for < t\ < ti < 1. 
We take $j>„ of the form (|3.18fl with 6(p instead of tp {6 > 0), and with smooth Xn approximating 
X = l[ti,t 2 ]- Using (|3.9() . (|3.7|) and ()3.6|) and letting n — ► oo we obtain 
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where 

i-T 



H T {9) = — / T T „ s (vXT{t-s)ve, T {;s)){x)dsii{dx), (3.20) 

t T jR d JO 

and vg^T is defined by (|3.8jl for *&t( x , t) = G^p( x )XT{t)- This vg t x also satisfies (|3.9|) . 

Unlike [S], where fourth moments were employed, we need moments of (Xt^) — XT(ti),(p) of 
arbitrary high order. By (|3.19|) we have 



jk 

E(X T (t 2 ) -X T ( tl )^) k = (-l) k — e H ^ 



A; = 1,2,.... (3.21) 



Using (ESJ) and (j3~2ll|) we have 

H T (0) = 0, 
^t(O) = 0, 

4 fc) (o) = (-i) fc ^/ d f P... fr^K^^.-W...)^) 

-T T JR d JO JO JO 

xxT(T-8fc)...XT(T-ai)dsi...ds fc /i(da;) ) fc > 2. (3.22) 
By (|3.21|) and (|3.22j) , tightness will be proved if we show that there exists 5 > such that 

|i4 fc) (0)| < C(k,<p){t 2 -h) kS for fc = 2,3,.... (3.23) 

The scheme described above is employed in the proofs of all results for 7 < d. The proof of each 
specific case, however, requires slightly different and non-trivial calculations; nevertheless, for brevity 
we will omit some proofs, concentrating on arguments which are either the most typical or the most 
involved. 

3.2 Proof of Theorem 2.1 

We will prove the theorem for fx(dx) = |x| _7 dx since in this case the formulas are slightly simpler. It 
will be obvious that the same type of argument applies for ^ of the form (2.1). It is easy to see that 
in this case the right hand side of 1)3. 14|) with $ given by (|3.16|) is of the form 

expj.Ki^y <p(x)dx^J y jf (u- s)- d/a s-^ /a x(s)x(u)dsd u y (3.24) 

Using (JXT2J), (ESS) and substituting u' = 1 - u/T, s' = 1 - s/T we obtain 

I ( T ) = p2 /_, / / i^PTs{x - y)f{y)p T (u-s)(y ~ z)f( z )x(s)x(u)\x\~' y dydzdsdudx. 



t JR Jo Jo 

We apply the self similarity of the a-stable density, i.e. 

Pat (x) = a-^ptixa- 1 ^), (3.25) 
substitute x' = xT~ l / a and use (2.2), then 

J(T) = I g T (x)\x\^dx, (3.26) 

.Iwd 
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where 

-1 l-U 



JO 



9t(x)= / / x(s)x{u)Ps{x-yT- 1 ' a ) Pu ^ s ({y-z)T- 1 l a )^(y)^{z)dydzdsdu. (3.27) 



By self similarity again the integrand in (|3.27|) is bounded by Cs d / a {u — s) d / a (p(y)tp(z), which is 
integrable since d < a. Hence, by the dominated convergence theorem we obtain 

1 i-u 

x(s)x(u)p s {x)pu-s(0)tp{y)<p{z)dydzdsdu 



lim g T (x) = 9oo (x)= / / 

r ^°° JO JO JR< 



= Pl(0) C rx{s)x{u)s- d ' a {u-s)- d l a p 1 {xs- 1 / a )dsdu( I <p(y)dy) . (3.28) 

Jo Jo \JR d J 

It is well known that 

PI{X) ~ l + |x|^ ' (3 - 29) 
hence, by (|3.25f) we easily deduce that for d < a, 

Jo Mx)ds ~ TTW^- 

This and an obvious estimate 



(3.30) 



1 + \ w \ d + a 

< C2 , T L„ (3.31) 



l + \x- w\ d+a ~ 1 + \x\ d+a 
(for w = yT~ l / a ) imply that 

9t{*) < C \ + l\ d+a - (3-32) 

By RFM . (EOfl) . (13321 and taking into account (l3~2l1) . we obtain (ETTil . 
This completes the proof of 1)3.31) . 

We proceed to the proof of (|3.15|) . We use ()3.13j) . ()3.10|) . ()3.18[) and boundedness of x to get 



C 



T 



H{T) < -^3 / / / / T T ^ s ((pT s _ u (ipT u ^ v <p))(x)dvduds\x\ 7 dx. (3.33) 



T JR JO JO JO 

„/ _ 1- 



We substitute v' = then v! = ^y^, then s' = -jA, and we increase the time intervals to [0,1], 
obtaining 



II(T) < f I f 

Ft J^- d J$t 3d Jo 



PTs(x - y)ds(p(y) p Tu (y - z)ducp(z) p Tv (z - w)dv 
Jo Jo 

X(p(w)dydzdw\x\~ 1 dx. (3.34) 

Denote 



f(x)= [ p s (x)ds (3.35) 

JO 



and 

£ T (x) =T d l a v{T l l a x). (3.36) 
Note that / is integrable, and by 1)3. 30 j) it is bounded, and 

(fT(x)dx = / ip(x)dx. (3.37) 
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Using ipPHJ) . substituting x' = xT~ l l a , y' = yT' 1 ^,^ = zT- l l a ,w' = wT~ l l a , we write (l3~3H) 

Tn3-(2d+ 7 )/a /• 

/I(T) < C 3 / a T (x)\x\-^dx, (3.38) 



T 

where 

ot(x) = f * ($T * (/ * (£t(/ * ^t))))(^)- 
The properties of / and y>r easily imply that 



sup / ax{x)dx < oo and sup sup o,t(x) < oo. 



Hence $TTB follows from (EOSl since 7 < d and T 3 -^ 2d+ ^/ a /F^ (see 

According to the general scheme, in order to prove tightness we show (|3.23|) . We substitute 
s'j = 1 — in 1)3.22(1 and we obtain 



|# (fc) (0)| (3.39) 
'-Zk / / ■•• / T T Sfc (^T( Sfc _ 1 - Sfe )(^--0--0(^)x(sfc)---X5(si)^i...ds fc — dx. 

-^V JR d J -'«■ 1^1 



= it!. 



We need the following estimate: 

Jr / V^(^M dM ^ CT^' 2 <*{t 2 - (3 40) 

(recall that x = l[ti,t 2 ])- By (|3-25() and boundedness of pi we have 

I' r r(w _ a)¥ >(y) X («)d« < CT^/ 2a J ^{z)dz j\u - s)- d ' a X {u)du. (3.41) 
Hence (|3~4"0|) follows. We iterate (|3~l0|) fc-1 times in (j3~39|) . estimate (T(T-«/ 2a ) fe - 2 by 1, arriving at 
|# (fc) (0)| < Cfa-t^-W-WTfr-W^^- f f f PTsk (x-yMy)dyd Sk \xrdx 

jR d Jo JR d 

= C(t 2 - /" /*^ T ( X )| X |-Tdx, 

where we have used (|2.2|l . self-similarity and the usual substitutions x' = xT -1 /", y' = yT~ 1//Q , where 
/ and are defined by (|3.35|) and (|3.36J) . Hence we obtain (|3.23() by the properties of / and </>y. 
The proof of the theorem is complete. □ 

3.3 Some properties of the a-stable semigroup in the critical case d = a. 

We will need the following facts, valid for d = a, tp € 5(R d ), cp > 0: 



1 



T 



sup sup - — — / T u tp(x)du < 00, (3.42) 

T>2 xm d logT J 

lim — ?— f TMx)du = Pl (0) [ ^{y)dy, (3.43) 

t^oo log t y y R d 

lim — — / / T u (p(x)du\x\~ d dx = 0. (3.44) 

T^oo logT J\ x \d >T Jo 



11 



These properties are perhaps known but we have not been able to find references for them, so we show 
briefly how to derive them. 

To prove (j3.42j) and ((3.43)1 . it is clear that it suffices to consider jT . We use self-similarity, then 
make a substitution which turns out to be particularly useful in the critical cases and will be applied 
several times. Namely, we put 

u = — , (3.45) 

obtaining 

-i- / T u cp(x)du = [ [ Pl ((x-y)T' u / d )ip(y)dydu. 

Hence 1)3.42(1 and ((3.43)1 follow immediately. 

To prove ((3.44(1 we again replace J Q T by and make the substitution ((3.45)1 . We then have 

T 

T u (p(x)du\x\~ d dx 



\x\ d dydudx, 



logT J\ x \d >T J 1 

= I C I P i{{x-y)T- u ' d My)V 
J\x\ d >TJo Jm d 

^ c ! f I i z i Lt-zS 1 + \y\ 2d T- 2u My)\x\~ d dydudx, 

J\x\*>tJo Ju d 1 + \x\ M T 2u 

where the last estimate follows from ((3.29)1 and (|3.31() (recall that d = a). As tp G 5(M rf ), this 
expression is estimated by 

d / \x\- 3d dx CT 2u du< ° 2 



\x\ d >T Jo lo g T 



by calculus. This proves ((3.44)1 . 



3.4 Proof of Theorem 2.3 



We will present only an outline of the proof. 

Following the general scheme, and again taking for simplicity fi(dx) = |x| _7 (ix, we prove that (see 

fl3~m (i3~mfl3~TS|n 

lim I(T) = K 1 [ s^ /a x 2 (s)ds[ [ ip(x)dx). (3.46) 

T ^°° Jo \jR d J 

In 1(3.12)) we substitute u' = T — u, s' = 1 — ^ , use self-similarity and put x' = xT~ l l a s~ l l a . By ((2.4)1 
we obtain 

I{T) = T^l I' f ^ x(s)x{s + ^]s^ /a Vi(x-ys- l l a T- 1 l a ) (3.47) 

X{p(y)T u ip(y)\x\~ 1 dudydsdx. 

Using 

sup / p\(x + z)\x\~ y dx < oo, (3.48) 



Zg 

it is easy to see that 



lim I(T) = lim I'(T), 
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where 

l'{T) = -i- / f T / / T(1 5) xWxf« + ^V" 7/a Pi( fl! -f a " 1/a:r " 1/a ) 

x (p(y)T u (p(y) \x\ _1 dudydsdx. 
We use self-similarity again and make substitution 1)3.45(1 . Then 

i\t) = [ [ T [ [ logT x^x^ + r"- 1 )^/^^-^- 1 /^- 1 /") 

JR d Jo JR d JO 

Xip(y) / pi((y - z)T~ u / a )ip(z)dzdudyds\x\~' y dx. 

JR d 

Now it is clear that the limit of I'(T) should have the form (|3,46|) , We omit details. 
We proceed to (|3.15|) . We use (|3.33() . substitute v' = u — v, then v! = s — u, then s' = and 
increase the time intervals appropriately, obtaining 

H(T) <C-^3 J T Ts (^p J T u(v T v Lpdv^du^j (x)ds\x\-^dx. (3.49) 

By (|3.42J) applied twice we have 

II(T) < Ci r ( lQ g r ) 2 f f T Ts ip(x)ds\x\-^dx. 
Ft J^ d Jo 

Hence (|3.15(l follows by self-similarity, substitution x' = x{Ts)~ l / a , ((3.48)1 and (|2.4[) . 

Tightness is proved similarly as in Theorem 2.1. Here are the main steps. Instead of ()3.38)) we 
show that ^ 

Jr j T T{u _ s)i p(y) X (u)du < C(t 2 - t x )\^l a \ 
we iterate this estimate k — 2 times in (|3.41)1 . and we obtain 

|# (fc) (0)| < Ci(t 2 -ti) (1 " 7/a)(fc ~ 2)/4 ^ /, (o). 

Using ()3.42|) it is not difficult to prove that 

H"(0) < C 2 (t x p ja - t{^ /a ) < C 2 (t 2 - h) 1 ^. 
Hence (EHE! follows. □ 

3.5 Proof of Theorem 2.4 

Again we give only a sketch of the proof. We make the same substitutions as in the beginning of the 
previous proof, and by ()2.6j) we obtain 

I{T) = Itl f T(1 S \{s)x(s + ^]s-^ a Pl {x-ys- l / a T- l / a My)T uV {y)\x\-^dudydsdx 

jR d Jo JR d Jo V 1 J 

(cf. 1)3.471) ). Hence it is not hard to see that 

lim I(T) = [ Pl (x)\x\~^dx [ s~ l/a X 2 {s)ds [ <p(y)G<p(y)dy, 

T^oo J R d J J R d 
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since ' T u ipdu = Gp. This implies (|3.14|) . Next, by (|3.49|) (which is always valid), 

II(T) < C^- I [ T Ts (^G(<pG<p))(x)ds\x\-~< dx. 
J^ d Jo 

Hence ()3.15j) follows by the usual argument since Gp is bounded. 

Tightness can be proved in the same manner, even easier, as in Theorem 2.3. □ 



3.6 Comments on the proofs of Theorems 2.5-2.7 

The proof of Theorem 2.5, though by no means straightforward, is slightly simpler than the proof for 
the doubly critical case (Theorem 2.8), which will be given in detail. Therefore we omit it. Proposition 
2.6 is obtained immediately from ()3.1|) and the following estimate, valid for ip > 0: 

* = £0*.)^* < c>l - - ; + \ x]1 < », 

where we have used the fact that sup x (l + |x| rf_a )|G</?(x)| < 00 ^3] (Lemma 5.3). 

The proof of Theorem 2.7 is similar (but not identical) to the argument carried out for Theorem 
2.1. We omit it for brevity. 



3.7 Proof of Theorem 2.8 

We apply the general scheme. By (|3,12|) . (J3.16|) - (|3.18|) and the substitutions v! = s — u, then s' = T—s, 
we have 

I(T) = h{T) + I 2 (T) + J 3 (T) + / 4 (T), (3.50) 



where 



h(T) = ^Lm^tL I "W.WMx^)^^)^*^, (3.51) 
1 



T rT-s r T-l rT- 



h(T) = -j / / / -/ / )..., (3.52) 

*T Jl<\x\ d <T \Jo JO Jl Jl 

h(T) = i / [ T [ T ~ S ..., (3.53) 

*T J\x\ d >T JO Jo 

h(T) = 4/ f r~ S -.., (3.54) 

*T J\x\<lJ0 Jo 



where . . . denotes the same integrand as in Ii(T). 
We will show that 



lim IAT) = K lX 2 (0) ( [ if(x)dx) , (3.55) 



and the remaining integrals converge to 0. 
By (t2~T2l 

h(T) < -7-. — =^7 / I / / . . . duds + / / ... duds ) dx. 

(logTy y N d< r \J Q Jo 



Jo 



Using (|3.42|1 we obtain 



C f 1 C 

h(T) < t- — ^—75- / r-rrcte < — 4^ as T -> 00. (3.56) 

2V ' ~ (logT)2 J ]x]d < T l + \x\* - logT 1 ; 
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The fact that 



lim I 3 (T) = (3.57) 

T — >oo 



follows immediately from (|3.42|) and ()3.44j) . and 

lim h{T) = (3.58) 

T — >oo 

is also a consequence of (|3,42j) . 

By (I33T1) . ()2~T2l and (ET23I) we have 

h(T) = --^—.t [ I 1 I* S S - 1 Pl ((x-y) S - 1 / d My)u- l p l ((y-z)u- 1 / d ) 



(logT) 3 Ji<\ x \d< T J R 2d J i 



X <f{z)x [ 7£ J X ( ^ + ^ J 1 \< d dudsdydzdx. 



We make the substitution (|3,45|) for both u and s, obtaining 

log(T-l) WT-Tf) 

Ji(T) = 7—^/ / / losT / losT pi((x- y )r- s /V((y-^- u/d ) 

x V 9(y)^(z)x(T s - 1 ) X (r s - 1 + T u ~ 1 )—^——^dudsdydzdx. 

1 + |x| a 

In the integral f dx we pass to polar coordinates (r,w) (r = \x\) and then substitute r' = r d . 
We have 

h(T) 

logTdJi Js^Jr^Jo Jo [U ' 1 1°. logT i 

x i p(y)f(z)x(T s ~ 1 )x{T s ' 1 + T u - 1 )-^—dudsdydza d ^ 1 {dw)dr, (3.59) 

1 + r 

where cx^-i is the Lebesgue measure on the unit sphere Sd-i in M. d . Again, we use ()3.45j) putting 
r' = log rj logT, then it is easy to see that the integrand converges to 

l[0,i] (s)l[o,i] (n)l [0 , s] (r)pj(0)cp(y)<p(z) X 2 (0), 

and is bounded by Cp\ (0)ip(y)(p(z). Hence (1335)1 follows. By 1)330)) and l)335)) - l)33%) ) we obtain (13141 . 

Now we pass to the proof of ()3.15|) . By (|3.13|) and ()3.1U|) for $ of the form (|3.16|) , after obvious 
substitutions we have 

T 

dx. 



Using (|3.42|) twice we get 

C f f T 1 
U{T) < -, rTTTTT / / T s <f(x)ds -—;dx 

v ' - (i og T)5/2y Md y > i + \ x \d 

= A + B, 

where 

a=-. ° f ...< c ' 



x\ d 



(l0gT)5/2 J\ x \ i>T '~- (logT)3/2 
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by (JS3U), and 



C f C f 1 

B = - -— / . . . < 2 -— - / r^jdx, (3.60) 

(logT)V2 J ]x]d < T " (logT)3/2y Nd < T l + |x| rf 1 ; 

< c * 



(log T) Va- 
in the first estimate in (3.60) we have used ()3.42j) once more. Hence ()3.15(l follows. 

Passing to the proof of tightness, first observe that the method employed in the proof of (|3,4|) 
can be also used to obtain convergence of finite dimensional distributions of Xt- This fact has been 
already used in jSj; here we repeat briefly the argument. For cp±, tp2, . . . ,<pk G 5(IR a! ), all ipj > 0, and 
< ti < . . . < tk < 1, it is easy to see that -Eexp{— Ylj=i(-^-T(tj), Pj)} has the form (|3.7|) with 

k 
3=1 

and the corresponding vt given by (|3.8|) , 

Approximating ^ by smooth functions we obtain that ()3.9j) holds, and then we argue as before. 

In particular Xt{s) converges in law. Therefore, to prove tightness of X T in C([e, 1],<S'(R )) it 
suffices to show that (X T - X T (e), ip) => in C([e, 1]), for any p G S(R d ), (p>0. 

Denote 

1 f tT 

wr(t) = -=- / {N a ,<p)ds, t > e. 

*T JeT 

By (1.2) it is clear that it is enough to show that wt and Ewt converge to in law in C([e, 1]). Both 
processes are increasing, so it suffices to prove that 



lim Ew T (I) = 0. 

T — >oo 



By (j3.1|) and substitution x' = xs 1 l a , 



Ew T (l) = — / pi{x-ys 1/a )p{y) , dxdyds 

F T JeT JR 2d 1 + M « 

= Ji(T) + J 2 (T), 



where 



b T JeT J\x\<\ 



?T JeT J\x\>\ 



We have 



J 2 (T) < £ C \- x ds = C 

TT JeT 

Ji(T) < / / ri ^dxds 



log(l/g) 
F T 



Ft JeT J\x\<l 1 + \ x \ s 
°* ^ 1 -drds<^m^0. 



Ft JeT Jo l + rs ~ (logT) 1 ^ 

□ 
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3.8 Systems with more general intensity measures \i 

In this section we consider a measure \x of the form ()2.13() . We sketch the proof of Proposition 2.9 and 
we discuss Example 2.10. 

Proof of Proposition 2.9 We concentrate on the case 7 < d < a. The other cases will be mentioned 
later. 

First notice that it suffices to assume that v = in ()2.13j) . since it is easy to see that with our 
normalization the terms corresponding to v will vanish in the limit. We repeat the steps of the proof 
of Theorem 2.1. Observe that boundedness of h implies that (j3,33|) also holds in the present case, 
hence (|3.15|) is obtained in the same way as before. Also, the tightness is proved without any change. 

It remains to show (j3.14|) . Instead of l)3.26j) we have 

I(T) = I grW ^^^ KT 1 '"*)**, (3.61) 

where gx is defined by ()3.27j) . We write 

J(T) = h(T) + J 2 (T), (3.62) 

where 

h(T)= f g 00 (x)-\-h(T 1 / a x)dx, (3.63) 

with given by (|3.28|) . and 

7 2 (T) = f ( g T (x) r7/ ° , - I h{T 1 / a x)dx. 

(l3~32|) implies that lim I 2 {T) = 0. 

T — *oo 

Note that from assumption ()2.14|) it follows that 

lim / a(x)h{Rx)dx = C a(x)dx (3.64) 

R^oo J R d J R d 

for a(x) = lui< r , where C is the limit (|2.14|) . Hence, it is easy to see that (|3.64|) also holds for 
any symmetric integrable function a. The function g o{x)\x\~' 1 is obviously symmetric and integrable 
(7 < d and ()3.32j0 . Therefore (|3,64j) implies (|3.14j) . This completes the proof in the case 7 < d < a. 

In the remaining cases for 7 < d, tightness and ()3.15j) follow immediately from the cooresponding 
proofs for [i of the form (|2.1|) . Also, to obtain (|3.14j) we repeat the same steps, obtaining I(T) in an 
analogous form as in H3.61j) . And then we apply ()3.64|) . □ 

For d = 7 the method described above cannot be applied. For example, in the case d = 7 = 1 < a 
the function goa{x)\x\~ d is not integrable. To prove (|3.14|) we would need existence of the limit 

1 f T l / a 1 

lim — / Qoo(x) ■ —r-hixT a)dx. 

T-oologT 7 N <i 'l + lxT 1 /^ V ; 

This limit is easy to obtain for h = 1 , but it is not clear how to formulate an elegant condition assuring 
its existence in a more general case. 
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The case 7 = d = a is even more complicated because in (|3.59|) we would have /i(wr 1//a ) under the 
integrals. 

Proof of non-existence of the limit in Example 2.10 It is obvious that the only nontrivial 
normalization is that given by (|2.2j) . since h is bounded and separated from 0. Analogously as in the 
previous proof, convergence of Xt is equivalent to convergence of I\{T) defined by (jUlESl). We will 
show that I\(T) does not converge. Let T n = n na+a /' 2 ,n = 2, 3, . . .. On the set {x : ^= < |x| < y/n} 
we have 



h(T^ a x) = u{n) = 
It is clear that lim {I\(T n ) — I[(n)) = 0, where 



1 if n is even, 

2 if n is odd. 



I' 1 (n) = u(n) / 9oo(x)\x\ 7 <ix, 



and obviously I[(n) does not converge. 
3.9 The finite measure case 



Proof of Lemma 2.12 We improve slightly the proof of Lemma 2 of Bingham [2]. Let 



It is easy to see using self-similarity that 



poo 

]im9 a ' / / e~ ds 



a' = 1 - -. (3.65) 
a 



(p(y)p s (x — y)dyds = K / tp(x)dx (3.66) 

Jr 



for any tp G <S(M), but in general the convergence is not uniform in x € R if (f is not compactly 
supported. Therefore condition (A) of Darling-Kac is not satisfied, so, unlike Bingham, we cannot 
apply directly their theorem. We prove that 

Z T j>KL\ (3.67) 

( j> denotes convergence of finite dimensional distibutions), where L is a continuous increasing process 
whose inverse is an a'-stable subordinator. On the other hand, it is known (see [Q, Prop. 4, Ch. V; 
see also |15| ) that such L is the local time process at of Q. 

It is clear that in order to prove ()3.67j) it suffices to show that 

{{Z(tx),(px),...,{Z(t k ),(p k ))=¥K[L(ti) [ i Pl (x)dx,...,L(t k ) [ Mx)dx) (3.68) 



for any ti,...,t k e [0, l],(pi,. £ S(R),(p 1} . . . , (p k > 0, k = 1, 2, . . .. 

Fix <pi, . . . ,(p k as above ((fj ^ 0). Let Mt be the measure on R + such that 



A- 



M T ([o, tl ] x ... x [o,t k ]) = F g n ( ^y - (3 - 69) 
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(|3.68|) will be proved if we show 

lim M T ([0, h] x . . . x [0, t k }) = EL{t\) . . . L(t k ) (3.70) 

T — >oo 

for all ti, . . . , t k £ [0, 1]. To this end, by Lemma 3 of [2] it suffices to prove that 
lim / e- {eitl+ - +ektk) M T (dt 1 ,...,dt k ) 

= [(**(!) + • • • + e *{k))( *w + • • • + V)) • • • e Ak)] ~~ a ' ( 3 - 71 ) 

for all 6i, . . . , 6 k > 0, the sumation being over all permutations tt of {1, . . . , k}. For simplicity we will 
show (|3.71|) for k = 2. Without loss of generality we may assume that (\,Lfj) = 1, j = 1,2. 
We have 

"I l-CO [ 

J e- eh - et2 M T (dt 1 ,dt 2 ) = (Ji(T) + J 2 (T)) , (3.72) 



where 



-i roc rt2 

Ji(T) = -2 / / e-^-^^EMCTtMCrt^dt^, (3.73) 

?T JO JO 

J 2 (T) = / / ...dt 2 dh, (3.74) 



T JO JO 



where . . . denotes the same integrand as in (|3,73|) , 

By the Markov property, self-similarity and (|2.17j) . we have 



JO JO JRJM. 

x Pl {T~ l l a t- 1/a y) Pl (T- l / a (t 2 - hy^iz - y))dzdydhdt 2 , 



hence 



lim MT) = ~ 

r ^°° Jo Jo q 



oo pt 2 e -6\t\ e -6 2 t 2 



-dt 2 dt lP {(0) 



= K 2 {9 X + 6 2 )- a 9z a , 

where K = — r(^)r(l — — ). The limit of J 2 (T) is calculated identically, so we obtain (|3.71|) . 

For (p > 0, (|3.67|) implies that (Zt,*p) => KL(\,(p) in C([0, 1]), since (Zt,^>) is an increasing 
process. From this it follows immediately that {(Zt, <£>)}t is tight for any (p £ 5(IR rf ). Hence the proof 
of the lemma is complete by Mitoma's theorem |24| . □ 

Proof of Theorem 2.11 It is easy to see that Lemma 2.12 remains true with the same limit if £ is 
replaced by x + (. 
On the other hand, 

tT 



1 F 

(Y T (t),<p) = jr J V(x 3 +Ql)ds 



xj eiVo 
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where £ , £ 2 , • • • are independent copies of £, independent of Nq. Now the theorem follows from Lemma 
2.12 and the fact that A^o(M) has the same law as ^. □ 

Proof of Theorem 2.13 Let Zt be defined by (|2,19|) with Ft = logT. It suffices to prove that 

(Z T (l),<p) =*j»i(0)pi(A,¥>> (3.75) 

for ip > 0. Indeed, (|^.75|) implies that f^ 2 ^ ip(£ s )ds converges to in probability for any < t\ < 
t2- Now it is easy to see that 

((Z T (tx),(px), . . . , (Z T (t k ),tp k )) px{0){px(\,<px), • • -,Pi(\<Pk)) 
for < tx < ■ ■ ■ < tk- Hence, analogously as before we obtain that Zt =>- KpxX in C([e, l],5'(IR d )), 



and this easily implies (|2,20|) , 

Observe that (|3.75|) has exactly the form as in the Darling-Kac theorem but we cannot apply 
it directly, since 

j roc p 

e~ 9s <p(y)Ps(x -y)dyds 



— lOg JQ 

may not converge uniformly in x as 9 — > (so, condition (A) is not satisfied). Nevertheless, the proof 
of the Darling-Kac theorem can be repeated with some care in the present case yielding the desired 
result. □ 

Appendix 

Proposition For 7 > 0, 

E(N t , <p) -> as t-*oo, f G S(R d ). 

Proof It suffices to assume <p > 0. 

Let 7 > d. Using self-similarity of the standard a-stable process, we have from (3.1), 

E(N t ,cp)= f f p t {x - y)ip(y)dy 



i ' v ' 1 + \x[y 

dx 



t- d/a / / Px{t- 1/a {x-y)My)dy 



1 + tar 7 ' 



and the integrals converge to a finite limit as t —* 00. 

Let 7 < d. By the following version of Young's inequality |2J (Theorem 4.2), 

<c|l/yisiwi%, 



f(x)g{x - y)h{y)dxdy 
p,q,r > 1, 1/p + 1/q + l/r = 2, we have 



£7(JV tj v)<c( [ pAxfdx] ^ [ (p(x)dx ( [ ^^dx^) ' 

where p > 1, q = l,r > d/7 > 1. The integrals are finite. Again by self-similarity, and the change of 
variable x = t l l a y, 

([ p t {xfdx\ P = t-(l-Vp)d/° ( /" p^yfV 
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where 1 — 1/p > 0. 

In both cases the result is obtained. 

□ 
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